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Transverse spectra of both jets and hadrons obtained in high-energy pp and pp collisions at 
central rapidity exhibit power-law behavior of 1/plf at high pp. The power index n is 4-5 for jet 
production and is 6-10 for hadron production. Furthermore, the hadron spectra spanning over 
14 orders of magnitude down to the lowest pp region in pp collisions at LHC can be adequately 
described by a single nonextensive statistical mechanical distribution that is widely used in other 
branches of science. This suggests indirectly the possible dominance of the hard-scattering process 
over essentially the whole pp region at central rapidity in high energy pp and pp collisions. We 
show here direct evidences of such a dominance of the hard-scattering process by investigating the 
power indices of UA1 and ATLAS jet spectra over an extended pp region and the two-particle 
correlation data of the STAR and PHENIX Collaborations in high-energy pp and pp collisions at 
central rapidity. We then study how the showering of the hard-scattering product partons alters 
the power index of the hadron spectra and leads to a hadron distribution that may be cast into a 
single-particle nonextensive statistical mechanical distribution. Because of such a connection, the 
nonextensive statistical mechanical distribution may be considered as a lowest-order approximation 
of the hard-scattering of partons followed by the subsequent process of parton showering that turns 
the jets into hadrons, in high energy pp and pp collisions. 
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I. INTRODUCTION 


Transverse momentum distribution of jets and hadrons 
provide useful information on the collision mechanisms 
and their subsequent dynamics. The transverse spec¬ 
tra of jets in high-energy pp and pp experiments at high 
Pt and central rapidity exhibit a power-law behavior of 
1 /pj, with the power index n ~ 4 - 5, which indicates 
that jets are scattered partons produced in relativistic 
hard-scattering processes dHl3|. On the other hand, the 
power index for hadron spectra are in the range of 6 to 
10, slightly greater than those for jets [ 3 . 191-0 . revealing 
that hadrons are showering products from jets, and the 
hadron spectra are modified from the jet spectra but re¬ 
taining the basic pow er-law structure of the jet spectra. 
It was found ll|-[l5| recently that the hadron spectra 
spanning over 14 decades of magnitude from the lowest 
to the highest pp at central rapidity can be adequately 
described by a single nonextensive statistical mechani¬ 
cal distribution that is widely used in other branches of 
sciences Em, 


F (pp) = A 1 




a) 


Such a distribution with g = l + l/nis phenomenologi¬ 
cally equivalent to the quasi-power law interpolating for¬ 
mula introduced by Hagedorn fist and others 0 

F( Pt ) = a( 1 + ^f) , (2) 

V PoJ 

for relativistic hard scattering. Both Eqs. m and © 
have been widely used in the phenomenological analysis 
of multiparticle productions, cf., for example, 0-0 and 
references therein. 

It is of interest to know why such a nonextensive statis¬ 
tical mechanical distribution © may be a useful concept 
for hadron production. It may also be useful to contem¬ 
plate its possible implications. The shape of the spec¬ 
trum reflects the complexity, or conversely the simplic¬ 
ity, of the underlying production mechanisms. If there 
are additional significant contributions from other mech¬ 
anisms, the specification of the spectrum will require de¬ 
grees of freedom additional to those of the relativistic 
hard-scattering model. The small number of apparent 
degrees of freedom of the spectrum over such a large 
Pt regiorQ suggests the possible dominance of the hard- 
scattering process over essentially the whole region of pp 
at central rapidity 00. 
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1 There are only three degrees of freedom in Eq. (1): A , q (or 
equivalently n), and T. Notice that three degrees of freedom are 
almost the minimum number to specify a spectrum. Our spec¬ 
trum is therefore very simple. It is interesting therefore that the 
counting of the degrees of freedom in our case can lead to the 
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The counting of the degrees of freedom provides merely 
an indirect evidence for the dominance of the hard- 
scattering process over the whole pt region. We would 
like to search for direct evidences for such a dominance in 
three different ways. The hard scattering process is char¬ 
acterized by the production of jets whose transverse spec¬ 
tra carries the signature of the power index of n ~ 4 - 5 
at central rapidity HHIm! . The relevant data come from 
well-defined jets with transverse momenta greater than 
many tens of GeV obtained in DO, ALICE, and CMS Col¬ 
laborations Ejl ET l . To seek direct supporting evidences 
of jet production by the hard-scattering process over the 
lower-pr region, we examine the experimental UA1 and 
ATLAS data which give the invariant cross sections for 
the production of jets from the low ~pt region (of a few 
GeV) to the high-p-r region (up to 150 GeV) [H, [H]. If 
the power index n of the UAl and ATLAS jet spectra at 
central rapidity is close to 4 - 5, it will constitute a direct 
evidence of the dominance of the hard-scattering pro¬ 
cess over the extended pr region, for pp and pp collisions 
at high energies. In such an analysis, we need to take 
into account important pr -dependencies of the structure 
function and the running coupling constant by refining 
the analytical formula of the hard-scattering integral. 

The hard-scattering process is characterized by the 
production of jets as angular clusters of hadrons. We can 
seek additional direct evidences for the dominance of the 
hard-scattering process by searching for hadron angular 
clusters on the near-side using the two-particle correla¬ 
tion data in high-energy pp collisions from STAR and 
PHENIX Collaborations 1311 - 11 Mi l . Two-particle angular 
correlation data are specified by the azimuthal angular 
difference A<j> and the pseudorapidity difference Arj of 
the two particles. If hadrons associated with a low- and 
high-pr trigger are correlated at (A</>, Ap) ~ 0 on the 
near-side, it will constitute an indication of the domi¬ 
nance of the hard-scattering process over essentially the 
whole pr region. 

Finally, the hard-scattering process is characterized by 
the production of two jets of particles. We can seek an 
additional direct evidence for the other partner jet by 
searching for angular clusters of associated hadrons on 
the away-side in two-particle correlation data from STAR 
and PHENIX Collaborations |dil Phil . A ridge of hadrons 
on the away-side at A<p ~ tt associated with a low -pr 
and high-pr trigger will indicate the production of the 
partner jet by the hard-scattering process over the whole 
Pt region. 

While our search has been stimulated by the simplicity 
of the hadron pt spectrum, it should be mentioned that 
the importance of the production of jets with pt of a 
few GeV (minijets) has already been well emphasized in 
the earlier work of 0 ,and the production of the low ~pt 


suggestion of possible hard-scattering dominance and the suc¬ 
cessful search for supporting direct evidences as we describe in 
the present work. 


jet (minijets) in the low ~pt region has been pointed in 
the work of [38l - [4ll ]. We are seeking here a synthesizing 
description linking these advances together into a sin¬ 
gle and simplifying observation on the dominance of the 
hard-scattering over the whole pt region, with a special 
emphasis on the production mechanism. Such a comple¬ 
mentary and synthesizing viewpoint may serve the useful 
purposes of helping guide our intuition and summarizing 
important features of the collision process. 

After examining the experimental evidences for the 
dominance of the relativistic hard-scattering process in 
the whole pt region, we would like to understand how jets 
turn into hadrons and in what way the jet spectra evolves 
to become the hadron spectra by parton showering. Our 
understanding may allow us to bridge the connection be¬ 
tween the hard-scattering process for jet production and 
its approximate representation by a nonextensive statis¬ 
tical mechanical distribution for hadron production. In 
consequence, the dominance of the hard-scattering pro¬ 
cess over the whole pt region may allow the nonexten¬ 
sive statistical mechanical distribution to describe the ob¬ 
served hadron transverse spectra spanning the whole pt 
region at central rapidity, in pp collisions at LHC. 

In this paper, we restrict our attention to the central 
rapidity region at p ~ 0 and organize the paper as fol¬ 
lows. In Section II, we review and refine the analytical 
results for the relativistic hard-scattering process. We 
use the analytical results to analyze the spectra for high- 
Pt jets in Section III. We note that jet spectra carry the 
signature of the hard-scattering process with a power in¬ 
dex n ~ 4 - 5 at central rapidity. In Section IV, we study 
the UAl and ATLAS data which extend from the low- 
Pt region of a few GeV to the high-pr region up to 150 
GeV. We find that the power index for jet production 
is approximately 4-5, supporting the dominance of the 
hard-scattering process over the extended pt region at 
central rapidity. In Section V, we seek additional evi¬ 
dences of the hard-scattering process from two-particle 
correlation data. In Section VI, we study the effects of 
parton showering on the evolution of the jet spectra to 
the hadron spectra. In Section VII, we examine the reg¬ 
ularization and further approximation of the relativistic 
hard-scattering integral to bring it to the form of the 
nonextensive statistical mechanics. In Section VIII, we 
analyze hadron spectra using the nonextensive statisti¬ 
cal mechanical distribution. We present our concluding 
summary and discussions in Section IX. 

II. APPROXIMATE HARD-SCATTERING 
INTEGRAL 

We would like to review and summarize the results of 
the hard-scattering integral obtained in our earlier works 
in 0 GM10 so that we can refine previous analyt¬ 
ical results. We consider the collision of A and B in 
the center-of-mass frame at an energy y/s with the de¬ 
tected particle c coming out at p ~ 0 in the reaction 
A + B —> c + X as a result of the relativistic hard- 
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scattering of partons a from A with parton b from B. 
Upon neglecting the intrinsic transverse momentum and 
rest masses, the differential cross section in the lowest- 
order parton-parton elastic collisions is given by 

E c d 3 a{AB —> cX) v 


dc 3 


dx a dx b G a/A (x a )G b / B (x b ) 
E c d 3 a(ab —$■ cX') 


ab 


dc 3 


( 3 ) 


where we use the notations in Ref. [12| with c the mo¬ 
mentum of the produced parton, x a and x b the forward 
light-cone variables of colliding partons a and b in A and 
B , respectively and da(ab^ cX') / dt the parton-parton 
invariant cross section. 

We are interested in the production of particle c at 
@CM = 90° for which analytical approximate results can 
be obtained. We integrate over dx a in Eq. m by using 
the delta-function constraint in the parton-parton invari¬ 
ant cross section, and we integrate over dx b by the saddle- 
point method to write 

[■ XaoG a / A (x a )][x b0 G b / B (xb )] = e^ Xb \ (4) 
We expand f{x b ) about its minimum at x b Q. We obtain 


af GbO 


) 9(x b o)i 


2tt 


~d 2 f{x b )/dxl\ Xi = Xb0 


For simplicity, we assume G a / A and G b / B to have the 
same form. At 0 C ~ 90° in the CM system, the minimum 
value of f(x b ) is located at 

Ct 


x b0 = x a o = 2x c , and x c = 


Vs' 


( 5 ) 


We get 

„ d 3 a(AB^cX) 

Ec - dV - 


B[x a0 G a / A (x a0 )][x b0 G b / B (x b0 )\ 
dcr(ab^cX') 


ab 


dt 


where 


B = 


2tt 


Tr(x b0 - c 2 T jx c s) V -d 2 f(x b )/dxl\ Xb =, 


•EbO 


For the case of x a G a / A (x a ) = A 0 (l — x a ) 9a , we find 


d 2 f(x b ) 


dx 2 


2ff(l - x c ) 


x c (l ir a o)(l x b Q ) 


( 6 ) 


( 7 ) 


( 8 ) 


and we obtaiifl 
_ d 3 cr(AB^cX) 

Ec - dV - 




(1 “ Xa,o) 9a,+ 2 (1 ~ X b o) 9b+ z 


ab V^9aV X c(l - Xc) 

da(ab^-cX') 

dt 


( 9 ) 


2 In Eq. (23) of the earlier work of [l3l , there was a typographical 
error in the quantity x^q in the denominator under the square 
root sign, which should be x c • 


The above analytical result differs from the previous re¬ 
sult of Eq. (9) of Ref. m , where the factor that appears 
in the above equation 

(1 - ZacOh 1 “ Xbo)*/\/0- ~ x c ) (10) 


was approximated to be unity. We wish to retain such a 
factor in order to obtain a more accurate determination of 
the power index, in cases where ct may be a substantial 
fraction of ^/s. 

If the basic process is gg —► gg, the cross section at 
9 C ~ 90° 0 is 


da(gg -> gg) 
dt 



( 11 ) 


If the basic process is qq' qq ', the cross section at 
0 C ~ 90° 0 is 


da(qq' —> qq 1 ) 47ra 2 5 

dt 9 16 


( 12 ) 


If the basic process is qq —> qq\ the cross section at 
9 C ~ 90° 0 is 


dcrjgg -A gg) _ 57ra|ll 
dt 4 Cj, 36 


(13) 


In all cases, the differential cross section varies as da(ab -A 
cX')/dt ~ a 2 /(c^) 2 . 


III. THE POWER INDEX IN JET 
PRODUCTION AT HIGH p T 


Our earlier investigation on the effects of multiple col¬ 
lisions indicates that without a centrality selection in 
minimum-biased events, the differential cross section for 
the production of partons at high-px will be dominated 
by the contribution from a single parton-parton scatter¬ 
ing that behaves as l/cf. 0,0-01 . It suffices to con¬ 
sider only the results of the single parton-parton collision 
as given in Eq. © which can be compared directly with 
the transverse differential cross sections for hadron jet 
and isolated photon production. 

From the results in the parton-parton cross sections in 
Eqs. 1)1111211311 . the approximate analytical formula for 
hard-scattering invariant cross section er; nv , for A + B —> 
c + X at g ~ 0, is 

d 3 cr(AB^cX) 
c dx? 

UC rj~0 

„ ot 2 {l-Xaoicr^+H 1 -x b0 (c T )) 9b+ ? /1/n 

CX - . = . : . (Id) 

CtV c t/VsV 1 ~ x c (c T ) 


We analyze the c T spectra by using a running coupling 
constant 


u s {Q{ct)) 


12?r 

27 ln(C + Q 2 /Aq CD ) ’ 


(15) 
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where Ajqcd is chosen to be 0.25 GeV to give a s (M|) = 
0.1184 [50|, and the constant C is chosen to be 10, both 
to give 0 ! s ((3~Aqcd) ~ 0.6 in hadron spectroscopy stud¬ 
ies [5l| and to regularize the coupling constant for small 
values of Q(c T ). We identify Q as c T and search for n by 
writing the invariant cross section for jet production as 

d 3 a(AB^cX) 

^"inv — j ~> 

dc r,~ o 

Aa 2 s (Q 2 (c T ))(l - XgQ (c T )) 9a+ ^ (1 - x b o(c T )) 9b+ s 
CjV 1 ^ x c (c T ) 

where the power index n for perturbative QCD has the 
value of 4.5. 

We identify the parton c with the produced jet and we 
define the jet transverse rapidity yp as the logarithm of 

ct/V~s, 

VT = ln • e,r = 7S' <17) 

then the results in Eq. m gives 

91n(Ti nv da s 2(g a + g b + l)e VT e VT 

-i - Tl —:-r I 1 o) 

dyr dyT 1 — 2e VT 2(1 — e VT ) 

and 

d 2 In cr inv _ d 2 a s 2(g a + g b + l)e VT e VT 

dy 2 . ~ dy 2 T (1 - 2e VT ) 2 2(1 - e^) 2 j 

Therefore in the (ln CT; n v)-(ln Et(ct)) plot in Fig. [ T ] the 
slope (<91n(7inv/d?/T) at small values of Et gives approx¬ 
imately the power index n and the second derivative of 
lncr inv with respect to In .Sr at large values of Ep gives 
approximately the power index g a +gb of the structure 
function. The exponential index g a = gb for the struc¬ 
ture function of a gluon varies from 6 to 10 in different 
structure functions [1, ]52j . We shall take g a = 6 from [hj. 

With our refinement of the hard-scattering integral in 
Eq. (9), our analytical invariant cross section of Eq. (fl6l) 
differs from that in our earlier work in [l2j] in the pres¬ 
ence of an extra energy- and p-r-dependent factor of Eq. 
(10) and a slightly different running coupling constant. 
We shall re-examine the power indices with Eq. m- 
We wish to obtain a more accurate determination of the 
power indices, in cases where ct may be a substantial 
fraction of the collision energy t/s. We also wish to use 
Eq. (il6l) to calibrate the signature of the power indices 
for jets at high p T , where jets can be better isolated, to 
test in the next section the power indices for jets extend¬ 
ing to the lower pt region, where jets are more numerous 
and harder to isolate. 

Using Eq. (USD, we find that the da/dgETdErlri^o data 
from the DO Collaboration |29| for hadron jet production 
within |r 7 |<0.5 can be fitted with n=4.39 for pp colli¬ 
sions at ^=1.8 TeV, and with n=4.47 for pp collisions 
at !/i=0.630 TeV, as shown in Fig. [Q 

In another comparison with the ALICE data for jet 
production in pp collisions at y/s = 2.76 TeV at the LHC 



FIG. 1. (Color online) Comparison of the relativistic hard- 
scattering model results for jet production, Eq. (1161) (solid 
curves), with experimental da/dyETdEr\ri~o data from the 
DO Collaboration [2§], for hadron jet production within 
|?7|<0.5, in pp collision at (a) Vs=L80 TeV, and (b) yjs=0.63 
TeV. 


within 1 77 1 < 0.5 [30| in Fig. 2, the power index is n=4.78 
for R = 0.2, and is n=4.98 for R = 0.4 (Table I). In Fig. 
3, the power index is n=5.39, for CMS jet differential 
cross section in pp collisions at -Js = 7 TeV at the LHC 
within |?j| < 0.5 and R = 0.5 31 . 



FIG. 2. (Color online) Comparison of the relativistic hard- 
scattering model results for jet production, Eq. (1161) (solid 
curves), with experimental da/dgErdET^^o data from the 
ALICE Collaboration (3(j, for jet production within |p|<0.5, 
in pp collision at 2.76 TeV for i?=0.4, and R= 0.2. 

The power indices extracted from the hadron jet spec¬ 
tra in DO [29|, (30|, and CMS (3lij are listed in Table I. 
The extracted DO power indices are smaller than those 
extracted previously in [l 2 j by 0.2 units, as the high¬ 
est transverse momenta are substantial fractions of the 
collision energy. In the other cases, the change of the 
power indices from our earlier work in [l 2 j] are small 
as their highest transverse momenta are substantially 
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smaller than the collision energies. 



FIG. 3. Comparison of the relativistic hard-scattering model 
results for jet production, Eq. (1161) (solid curves), with exper¬ 
imental da/dyEpdEp \ v ~o data from the CMS Collaboration 
[Uj, for jet production within |?7|<0.5, in pp collision at 7 
TeV. 


mechanism, extending the dominance of the relativistic 
hard-scattering process to the lower -pp region. 


IV. JET PRODUCTION IN AN EXTENDED 
REGION FROM LOW TO HIGH p T 

The analysis in the last section was carried out for jets 
with a transverse momentum greater than 19 GeV. It is of 
interest to find out whether the perturbative QCD power 
index remains a useful concept when we include also the 
production of jet-like energy clusters (mini-jets) at lower 
transverse momenta. In order to apply the power-law 
(El) to the whole range of ct , we need to regularize it by 
the replacement^. 


1 1 

Cj, 1 + rrij,/m^o 

or alternatively as 

1 1 

ct 1 + mp/mpo 


( 20 ) 


( 21 ) 


With the jet spectra at high pr from the DO, ALICE 
and CMS Collaborations, we find that the signature for 
jet production is a power index in the range of 4.5 to 
5.4, with a small variation that depend on yfs and R as 
shown in Table I. While these power indices are close to 
the lowest-order pQCD prediction of 4.5, there appears 
to be a consistent tendency for n to increase slightly as 
yfs and R increases. Such an increase may arise from 
higher-order pQCD effects. We can envisage the physical 
picture that as the jet evolves by parton showering, the 
number of generations of parton branching will increase 
with a greater collision energy yfs or a greater opening 
angle R. A greater yfs or a larger R value corresponds 
to a later stage of the evolution of the parton showering 
and they will lead naturally to a slightly greater value of 
the power index n. 


TABLE I. The power index n for the jet spectra in pp and 
pp collisions. Here, R is the jet cone angular radius used in 
the jet search algorithm. 




pp Region 

R 

V 

n 

D0|29J 

pp 1.80TeV 

64<pT<461GeV 

0.7 

^ROT 

4.39 

DO [29] 

pp 0.63Te"V 

22<p T <136GeV 

0.7 

M<0.7 

4.47 

ALICE[30 

pp 2.76Te'V 

22<p T <115CeV 

0.2 

H<0.5 

4.78 

ALICE[30 

pp 2.76TeV 

22<p r <115GeV 

0.4 

M<0.5 

4.98 

CMS [31] 

pp 7TeV 

19<pt < 1064GeV 

0.5 

p <0.5 

5.39 


The signature of the power indices for the production 
of jets at high pp can be used to identify the nature 
of the jet production process at low pp. If the power 
indices in the production in the lower-p^ region are sim¬ 
ilar, then the jets in the lower-pp region and the jets 
in the high-py region have the same spectral shape and 
can be considered to originate from the same production 


The quantity mp o measures the average transverse mass 
of the detected jet in the hard-scattering process. Upon 
choosing the regularization El), the differential cross sec¬ 
tion d 3 a(AB -A pX)/dydpp in (TTUll is then regularized 
as 

d 3 a(AB -A pX) 
dydp T 

^ ^(Q 2 (cT))(l-^ao(cT)) 9a + 1/2 (l^(Ebo(cT)) gt + 1/2 ^ 
[1 + rrip(cp) / yj 1 - Xc(ct ) 

We fit the inclusive UAl jet cross sections at rj ~ 0 
[32f as a function of the jet pp for pp collisions with the 
above equation, and we find that n=4.47 and toto=0-267 
GeV for pp collisions at ^=564 GeV, and n=4.73 and 
mp o=0.363 GeV for pp collisions at ^=630 GeV. 

The ATLAS pp spectra for pp collisions at 7 TeV also 
extend to the region of a few GeV. It is of interest to find 
out what are the power indices for these collisions. We 
show in Fig. [5] the comparison of the results of Eq. Cdl) 
with the ATLAS data at 77 ~ 0 (33|. We find that n=5.03 
for R = 0.4 and n = 5.29 for i?=0.6. Because the data 
start with pp of a few GeV, the fits and the extracted 
value of n are insensitive to the variation of the TOto 


3 So far, the only rationale behind this is that, in the QCD ap¬ 
proach, large cp partons probe small distances (with small cross 
sections). With diminishing of cp , these distances become larger 
(and cross sections are increasing) and, eventually, they start to 
be of the order of the nucleon size (actually it happens around 
cp — cp 0 ~ l / r nuc ieon or m T — mpo )• At that point the cross 
section should stop rising, i.e., it should not depend anymore on 
the further decreasing of transverse momentum cp. The scale 
parameter mp 0 can then be identified with mp 0 here. Similar 
idea was employed when proposing Eq. G3- 
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FIG. 4. (Color online) Comparison of the relativistic hard- 
scattering model results for jet production, Eq. (1221) (solid 
curves), with experimental da/dr)ppdpp data from the UA1 
Collaboration [33], for jet production within |»y|<1.5, in pp 
collision at (a) -^5=0.546 TeV, and (b) -^5=0.63 TeV. 


TABLE II. The power index n extracted from jet production 
in pp and pp collisions in the extended pp region from a few 
GeV to the high-py region. 



Vs 

Pt Region 

R 

V 

n 

UA1|32| 

pp 0.564TeV 

5.5<p T <150GeV 

0.75 

M<l-5 

4.47 

UAl [32] 

pp 0.63TeV 

5.5<p T <150GeV 

0.75 

|??<1.5 

4.73 

ATLAS [33] 

pp 7TeV 

4.5<pT<95GeV 

0.4 

M<0.5 

5.03 

ATLAS [33] 

pp 7TeV 

4.5<pT<95GeV 

0.6 

M<0.5 

5.29 


few GeV (minijets) has already been emphasized in the 
earlier work of [ 3 ] • 

By comparing the power indices obtained in Table I 
for DO, ALICE, and CMS for jets at high pp with those 
for UAl and ATLAS for jets in the lower-pp region in 
Table II, we note that these power indices are very sim¬ 
ilar. The corresponding power index values are nearly 
the same, and the changes of the power index with re¬ 
spect to y/s and R are nearly the same. They can be 
considered to originate from the same relativistic hard- 
scattering mechanism, indicating the dominance of the 
hard-scattering process over the extended pp region from 
a few GeV to about 100 GeV. 



FIG. 5. Comparison of the relativistic hard-scattering model 
results for jet production, Eq. (1221) (solid curves), with experi¬ 
mental da/dr/prdpT\ri~o data from the ATLAS Collaboration 
[33| , for jet production within |p|<0.5, in pp collision at y/s=7 
TeV, with (a) 77=0.4 and (b) R= 0.6. 


values so that there is an ambiguity in the product of the 
normalization and toto in the analysis. The fits in Fig. 
J5]) have been obtained with mpo = 1 GeV. The value of 
n is is related to the slope of the curves in Fig. El 

We list in Table II the power indices extracted from 
UAl and ATLAS for the extended pp region from a few 
GeV to the high-pp region. It should be mentioned that 
the importance of the production of jets with pp of a 


V. ADDITIONAL EVIDENCES OF JET 
PRODUCTION FROM TWO-PARTICLE 
CORRELATION DATA 


In addition to the spectral shape, we seek addi¬ 
tional evidences of jet production in the low -pp re¬ 
gion from experimental two-particle correlations mea¬ 
surements, which consist of correlations on the near-side 
and the away-side. We shall first examine near-side cor¬ 
relations . 

The experimental distribution of near-side particles as¬ 
sociated with a trigger particle of momentum p^' 6 in pp 
collisions can be described well by [13, l43j 


dm 

_Jfi_ =N- 

ppdpp dAi] dAcf> J0t 


exp{ (m — sjm 2 +p^)/T- ]et } 


'2? xR 2 


T]et{jn + Tjet) 
o -[(A0) 2 + (A^) 2 ]/2ii 2 


(23) 


where by assumption of hadron-parton duality m can be 
taken as the pion mass m v , Nj et is the total number of 
near-side (charged) associated particles in a pp collision, 
and Tj e t is the jet inverse slope (“temperature”) param¬ 
eter of the “pp jet component”. We find that the param¬ 
eters TVjet and Tj e t vary linearly with p^ Ig of the trigger 
particle which we describe as 

Ajet = Ajeto + dpf Pp' S , (24) 


Tjet — T 


jetO 


dr Pp is - 


(25) 


We also find that the width parameter R depends slightly 
on pp which we can parametrize as 

R = R 0 - . (26) 


V m l + 


~i>p 
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Using this set of parameters and Eq. (l23l) . we fit the 



■10 1-10 1-10 1-10 1 


A(f) (rad) 


FIG. 6 . PHENIX azimuthal angular distribution of associ¬ 
ated particles per trigger in different p) rlg ® pj) ssoc combina¬ 
tions. The open circles are the associated particle yields per 
trigger, dN c h/ NtiigdAtp, in pp collisions (37j. The solid curves 
are the theoretical associated particle yields per trigger cal¬ 
culated with Eq. ((23j) . 


pp associated particle data obtained in PHENIX mea¬ 
surements for pp collisions at y/ snn = 200 GeV. The 
values of the parameters are given in Table III. As ex¬ 
tracted from Fig. 1 of (43j . the theoretical results of 
/N tlig dA(j> from Eq. C3l) are given as solid curves 
in Fig. 6 and the corresponding experimental data are 
represented by open circles. As one observes in Fig. 
6 , although the fit is not perfect, the set of parame¬ 
ters in Table III adequately describe the set of pp as¬ 
sociated particle data for 2 < p!)7 g < 10 GeV and 
for 0.4 < pf? soc < 4 GeV. As indicated in Table III, 
the parameters of Eqs. (l24l) and (l25l) are Aj e to = 0.15, 
d N = 0.1/GeV, Tj et o = 0.19 GeV, and d T = 0.06. It 
is interesting to note that the cone angle Rq for jets in 
the lower-px region is of the same order as those in the 
high -pp region. 

The presence of a well-defined cone of particles asso¬ 
ciated with a pr > 2-3 GeV triggers in Fig. 6 on the 
near-side and the non-vanishing extrapolation of the jet 
yield lVj e t to the case of a low-px trigger in Eq. (1!H1) 
provide an additional evidence of jet production in the 
px trigger > 2 GeV region in high energy pp collisions. 
Furthermore, even in minimum-px-biased events without 
a high-px trigger, a similar cone of associated correlated 
particles at (A</>, A 77 ) ~ 0 are present in t wo-p article cor¬ 
relation data, as shown in Fig. [7] |35|, f38M40l | . indicating 
the production of jet-likc structure on the near-side for 
low-px particles. 


TABLE III. Jet component parameters in Eq. (12311 obtained 
for the description of experimental near-side associated par¬ 
ticles with different p) rlg triggers in STAR [34| and PHENIX 
[37| Collaborations, in pp collisions at v /sjviv=200 GeV. 



STAR 

PHENIX 

Pr B 

4-6GeV 

2-3GeV 3-4GeV 4-5GeV 5-10GeV 

Ajet 

0.75 

0.15+0.10 (Pt 19 )/GeV 

Tjet 

0.55GeV 

0.19 GeV+0.06 (p" 19 ) 

Ro 

0.50 

m a 

1.1 GeV 



FIG. 7. (Color online) Minimum-pT-biased two-particle an¬ 
gular correlation, without a pr trigger selection, for charged 
hadrons produced in pp collisions at y7s=200 GeV (from Fig. 
1 of Ref. f3f| of the STAR Collaboration). Here, </>a is Ac/), 
the difference of the azimuthal angles of two detected charged 
hadrons, and t]a is At], the difference of their pseudorapidi¬ 
ties. 


In addition to the particles associated with trigger par¬ 
ticle on the near-side, there are particles associated with 
the trigger particle on the back-to-back, away-side at 
A()>~7r, in the form of a ridge along the Ap direction, both 
with high-px [dll [iibl and low-px triggers dM FiSl Fift , for 
pp collisions at yfs= 200 GeV. Here, the importance of 
the production the low-px jet (minijets) in the low-px re¬ 
gion has already been pointed out previously in the work 
of dBjdlJ. In Fig. 7, (taken from the STAR data in Fig. 
1 of [39|), we show the two-particle correlation data in a 
minimum-px-biased measurements which corresponds to 
the case with a low-px trigger. The two-particle correla¬ 
tion data in Fig. [7] indicate the presence of (i) a near-side 
particle cluster at (A cf>, Ap) ~ 0 (a mini-jet) and (ii) an 
away-side ridge of associated particles at A <f> ~ n. The 
A 4> ~ it (back-to-back) correlation in the shape of a ridge 
indicates that the two particles are parts of the partons 
from the two nucleons and they carry fractions of the lon¬ 
gitudinal momenta of their parents, leading to the ridge 
of Ap at A </> ~ 7 r. These two group of particles at A</> ~ 0 
and Acj) ~ 7 r can be interpreted as arising from the pair 

























































































of scattered partons in a relativistic hard scattering. 

The dominance of the hard scattering in the spectrum 
does not imply the absence of soft processes. It only stip¬ 
ulates that the soft process contribution is much smaller 
in comparison. In the lowest px region, one expects con¬ 
tributions from soft nonperturbative QCD physics that 
may involve the parton wave functions in a flux tube 
[53( | . the thermodynamics and the recombination of par- 
tonsjll, , or the fragmentation of a QCD string 

[57H59I ]. However, as the contributions from the hard- 
scattering processes increase with increasing collision en¬ 
ergies, the fraction of the contributions from soft pro¬ 
cesses becomes smaller in comparison with the contribu¬ 
tions from the hard-scattering processes, as pointed out 
earlier in 0 m. As a consequence, the contributions 
from the hard-scattering process can dominate the parti¬ 
cle production process in high-energy pp and pp collisions. 


VI. EFFECTS OF PARTON SHOWERING ON 
TRANSVERSE DIFFERENTIAL CROSS 
SECTION 


The last sections show the possible dominance of jet 
production at central rapidity in high-energy pp and pp 
collisions over essentially the whole px region. We would 
like to find out how the jets evolve to become hadrons 
and how the hadron spectra manifest themselves. 

In addition to the jet transverse spectra, experimental 
measurements also yield the hadron spectra without the 
reconstruction of jets. The hadron transverse spectra 
give a slightly greater power index, rihadron^6-10 0,S- 
Il4| . Previously, we outline how the increase in the power 
index n from jet production to hadron production may 
arise from the subsequent parton showering that turns 
jets into hadrons [l^l - We would like to describe here 
the evolution in more details. To distinguish between 
jets and its shower products, we shall use the symbol 
c to label a parent parton jet and its momentum and 
the symbol p to label a shower product hadron and its 
momentum. 

The evolution of the parton jet into hadrons by par- 
ton showering has been described well by many models 
[60[ . There are three different parton showering schemes: 
the PYTHIA [HI, We HERWIG [oil . and the ARIADNE 
[HI. The general picture is that the initial parton is char¬ 
acterized by a momentum and a virtuality which mea¬ 
sures the degree of the parton to be off-the-mass-shell. 
The parton is subject to initial-state and final-state ra¬ 
diations. After the hard scattering process, the parton 
possesses a high degree of virtuality Q^, which can be 
identified with the magnitude of the parton transverse 
momentum cx- The final-state radiation splits the parton 
into binary quanta as described by the following splitting 


DGLAP kernels [64|. 


P _H + * 2 

q ^ qg "31 - * ’ 

(27) 

9 z(l-z) ’ 

(28) 

P g ^ = ^[z 2 + (l-z)% 

(29) 


where z is the momentum fraction of one of the showered 
partons, and there is symmetry between z and 1 — z for 
symmetrical products in the second and third processes. 
After the showering splitting processes, there is always a 
leading parton with 


^leading iA ^non-leading- 


(30) 


For the study of the px hadron spectra as a result of 
the parton showering, it suffices to focus attention on the 
leading parton after each showering splitting because of 
the rapid fall-off of the transverse momentum distribu¬ 
tion as a function of increasing cx ■ As a consequence, we 
can envisages the approximation conservation of the lead¬ 
ing parton as the parton showering proceeds and as its 
momentum is degraded in each showering branching by 
the fraction (z)= (heading)- In the present study of high- 
px particles in the central rapidity region, the parton c 
is predominantly along the transverse direction, and the 
showering of the produced hadrons will also be along the 
transverse direction. 

A jet parton c which evolves by parton showering will 
go through many generations of showering. If we label 
the (average) momentum of the i-th generation parton 
by c~x , the showering can be represented as cx —> cP —> 

cP —> cP cP=px- Each branching will kine¬ 

matically degrade the momentum of the showering par- 
ton by a momentum fraction, {z)=cp rl pcP. At the end 
of the terminating A-th generation of showering, the jet 
hadronizes and the px of a produced hadron is related to 
the cx of the parent parton jet by 


Pt 

ct 



CT 


(31) 


It is easy to prove that if the generation number A and 
the fragmentation fraction z are independent of the jet 
ct, then the power law and the power index for the px 
distribution are unchanged IT3 . 

We note however that in addition to the kinematic de¬ 
crease of cx as described by m, the showering gener¬ 
ation number A is governed by an additional criterion 
on the virtuality. From the different parton showering 
schemes in the PYTHIA [HI , the HERWIG |62| , and the 
ARIADNE [( 33)1 , we can extract the general picture that 


the initial parton with a large initial virtuality de¬ 
creases its virtuality by showering until a limit of Q cutoff 
is reached. The downgrading of the virtuality will pro¬ 
ceed as Qj et =Q(°) QW Q( 2 ) -> Q( 3 ) ... 
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g(A) = g cut °ff j with 


n(*+i) qW 

(0 = ^TvrT” and = ( 0 X - 


Q« 


Q jet 


(32) 


The measure of virtuality has been defined in many 
different ways in different parton showering schemes. We 
can follow PYTHIA [60| as an example. We consider a 
parton branching of a —> be. The transverse momentum 
along the jet a direction is 


b^= z( 1 — z)a 2 — (1 — z)b 2 — zc 2 . 


(33) 


If a 2 = [QW] 2 =virtuality before parton branching, and 
b 2 = c 2 = 0, as is assumed by PYTHIA, then 

b 2 T = [Q( i+1 )] 2 = z( 1 - z)a 2 = z(l - z)[QM} 2 . (34) 

So, if we identify the transverse momentum bj, along the 
jet axis as the square of virtuality [Q^ +1 ^] 2 after parton 
branching, the quantity z(l—z) measures the degradation 
of the square of the virtuality in each QCD branching 
process, 




(35) 


Thus, the virtuality fraction of Eq. © is related to 
<z( 1 - z)) by 


<0 = V<*(i-*)>- 


(36) 


As z is less than 1, (£) < (z) which implies that on the 
average the virtuality fraction (£) in a parton branching 
is smaller than the momentum fraction (z). As a conse¬ 
quence, the virtuality of the leading parton is degraded 
faster than its momentum as the showering process pro¬ 
ceeds so that when the virtuality reaches the cutoff limit, 
the parton still retains a significant fraction of the initial 
jet momentum. 

The process of parton showering will be termi¬ 
nated when the virtuality Q^ reaches the cutoff value 
Q( x ) =Q cutoff . at which the parton becomes on-the-mass- 
shell and appears as a produced hadron. This occurs 
after A generations of parton showering. The generation 
number A is determined by 


A = In 


Q 


cutoff 


Q jet 


ln(£). 


(37) 


There is a one-to-one mapping of the initial virtuality 
Q jet with the initial jet transverse momentum ct of the 
evolving parton as Q jet (cT) (or conversely CT(Q jet ))- The 
cut-off virtuality Q cutoff maps into a transverse momen¬ 
tum CTo=CT{Q cutoB )- Because of such a mapping, the 
decrease in virtuality Q corresponds to a decrease of the 
corresponding mapped ct■ We can infer from Eq. © 
an approximate relation between ct and the number of 
generations, A, 


Q 


,cutoff (cTo)^\ /, , ( CTO 


Thus, the showering generation number A depends on 
the magnitude of the jet momentum ct- On the other 
hand, kinematically, the showering processes degrades 
the transverse momentum of the parton ct to that of 
the px of the produced hadron as given by Eq. ( 1311 ) , de¬ 
pending on the number of generations A. The magnitude 
of the transverse momentum pt of the produced hadron 
is related to the transverse momentum ct of the parent 
parton jet by 


PL = {z) * = {z) ^)/ ln <«>. 

ct 


(39) 


We can solve the above equation for pt as a function of 
ct and obtain 


and conversely 


where 


Pt 

f ct 

ct o 

\Cto 

— = 1 

( El 

CTO 

\Cto 




1/(1 -/0 


A* = 


Mf) 
ln(£>' 


(40) 


(41) 


(42) 


In practice p (or equivalently, the cut-off parameter 
Qcouoff Qr Ct ° ^ j s a p aram eter that can be tuned to fit the 
data. As a result of the virtuality degradation and virtu¬ 
ality cut-off, the hadron fragment transverse momentum 
Pt is related to the parton momentum ct nonlinearly by 
an exponent 1 — p. 

After the showering of the parent parton ct to the 
produced hadron pt, the hard-scattering cross section for 
the scattering in terms of hadron momentum pt becomes 


d 3 cr(AB pX) d 3 a(AB —> cX) dcT 


dydp T 


dydcT dp T 


(43) 


Upon substituting the non-linear relation (TUI) between 
the parent parton moment ct and the produced hadron 
Pt in Eq. m, we get 


dCT 

dp T 


1 


Pt 


1 — P \ CTO 


J2jj_ 

1-p. 


(44) 


Therefore under the parton showering from c to p , the 
hard-scattering invariant cross section <Ji nv (pt) for AB —> 
pX for hadron production becomes 


o-inv(pr) = E t 


d 3 cr(AB^pX) 


dp 3 


d 3 a(AB —y pX) 


y~ o 


dydp T 


y~ o 


a 2 s (Q 2 {c T ))( 1 - x a0 (c T )) 9 “ + 2 (1 - x b0 (c T )) 9f>+ 2 


where 


Pt V 1 - x c (c T ) 


, n — 2p . , ,1 

n = -. with n = 4 -)—. 

1-p 2 


(45) 


A = In 


MO- (38) 


(46) 
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Thus, the power index n for jet production can be sig¬ 
nificantly changed to n' for hadron production because 
the greater the value of the parent jet Ct, the greater the 
number of generations A to reach the produced hadron, 
and the greater is the kinematic energy degradation. By 
a proper tuning of p, the power index can be brought to 
agree with the observed power index n' in hadron pro¬ 
duction. The quantity p is related to n' and n by 


n' — n 



(47) 


For example, for /i=0.4 one gets n'=6.2 and for p = 0.6 
one gets n'=8.2. 


VII. REGULARIZATION AND FURTHER 
APPROXIMATION OF THE 
HARD-SCATTERING INTEGRAL 


In order to apply the power-law (H5l) to the whole range 
of pt for hadron production, we need to regularize it. 
Upon choosing the regularization m , the differential 
invariant cross section <7i n v (pt) for the production of a 
hadron with a transverse momentum p t becomes 


CTinv(PT) 


d 3 a(AB -A pX) 
dydp T 


Q«(Q 2 (ct))(1 ~ Xao(c T )) 9a+ Hl ~ x b o(c T )) 9b+1 2 ^ 

(1 + m T /m T o) n ' a / 1 - x c (c T ) 


In the above equation, the variable ct(pt) on the right- 
hand side refers to the transverse momentum of the par¬ 
ent jet ct before parton showering as given by Eq. (HD), 


qApt) _ f PT Y~ 2)/{n ~ 2) 

CTO \ CTO J 


(49) 


The quantities x a o, Xbo, and x c in Eqs. (l45l) are given by 
Eq. ©. 

We can simplify further the pt dependencies of the 
structure function in Eq. (1481) and the running coupling 
constant as additional power indices in such a way that 
will facilitate subsequent phenomenological comparison. 
We can cast the hard-scattering integral Eq. H51) for 
hadron production in the nonextensive statistical me¬ 
chanical distribution form 


d 3 a(AB -A pX) 
dydp T 


CTin v(Pr) 


A 

[1 + tot/toto]” ’ 


(50) 


where 


n = n' + riA, 


(51) 


and n' is the power index after taking into account the 
parton showering process, tia the power index from the 
structure function and the coupling constant. We con¬ 
sider the part of the pT-dependent factor in Eq. (l48l) 


f(pr) 


A s (c t (pt)){ 1 - 2c T {pT)/Vs) 9a+9b+1 
[1 ^ c t (pt)/V~A /2 


(52) 


that is a known function of pt- We wish to match it to 
a nonextensive statistical mechanical distribution with a 
power index n/\, 

KpA (1 + tot(pt)/tot 0 )" a 
We match the two functions at two points, pt i and pT 2 , 
f(PTi) = f{PTi), i = 1,2 (54) 

Then we get 

nA = _ ln/(pri)-ln/(pT2) _. (55) 

ln(l + tot(pt2)/toto)-1h(1 + m T (pTi)/m T o)' V 

As /(pt) is a known function of pt and y/s, «a can 
in principle be determined. The total power index n as 
given by (1511) is also a function of y/s. 

In reaching the above representation of Eq. (15UD for 
the invariant cross section for hadrons, we have approxi¬ 
mated the hard-scattering integral er; nv (pT) that may not 
be exactly in the form of A/[ 1 + mr/mTo]™ into such 
a form. It is easy then to see that the upon match¬ 
ing (Ti nv (pT) with A/\l + rriT/mTo] n according to some 
matching criteria, the hard-scattering integral erj nv (pt) 
will be in excess of A/[ 1 + mr/rn'ro] 11 in some region, 
and is in deficit in some other region. As a consequence, 
the ratio of the hard-scattering integral <7i nv ( pt) to the 
fitting A/[ 1 + mr/mTo]™ will oscillate as a function of 
Pt- This matching between the physical hard-scattering 
outcome that contains all physical effects with the ap¬ 
proximation of Eq. (1501) may be one of the origins of the 
oscillations of the experimental fits with the nonextensive 
distribution (as can be seen below in Fig. 8). 


VIII. SINGLE-PARTICLE NONEXTENSIVE 
DISTRIBUTION AS A LOWEST-ORDER 
APPROXIMATION OF THE 
HARD-SCATTERING INTEGRAL 

In the hard-scattering integral Eq. (l50l) for hadron in¬ 
variant cross section at central rapidity, if we identify 

1 T 

n -A- and toto -— nT, (56) 

q -1 9-1 

and consider produced particles to be relativistic so that 
rriT ~ Et ~ pt, then we will get the nonextensive distri¬ 
bution of Eq. (1) as the lowest-order approximation for 
the QCD-based hard-scattering integral. 

It is necessary to keep in mind that the outlines lead¬ 
ing to Eqs. (l48l) and (l50l) pertains only to average values, 
as the stochastic elements and distributions of various 
quantities have not been properly taken into account. 
The convergence of Eq. (15U1) and Eq. (1) can be con¬ 
sidered from a broader viewpoint of the reduction of 
a microscopic description to a single-particle statistical- 
mechanical description. From the microscopic perspec¬ 
tive, the hadron production in a pp collision is a very 
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complicated process, as evidenced by the complexity of 
the evolution dynamics in the evaluation of the pp spec¬ 
tra in explicit Monte Carlo programs, for example, in 
[6li - l63| . There are stochastic elements in the picking of 
the degree of inelasticity, in picking the colliding parton 
momenta from the parent nucleons, the scattering of the 
partons, the showering evolution of scattered partons, the 
hadronization of the fragmented partons. Some of these 
stochastic elements cannot be definitive and many differ¬ 
ent models have been put forth. In spite of all these com¬ 
plicated stochastic dynamics, the final result of Eq. (1501) 
of the single-particle distribution can be approximated 
to depend only on three degrees of freedom, after all is 
done, put together, and integrated. The simplification 
can be considered as a “no hair” reduction from the mi¬ 
croscopic description to nonextensive statistical mechan¬ 
ics in which all the complexities in the microscopic de¬ 
scription “disappear” and subsumed behind the stochas¬ 
tic processes and integrations. In line with statistical 
mechanics and in analogy with the Boltzmann-Gibbs dis¬ 
tribution, we can cast the hard-scattering integral in the 
nonextensive form in the lowest-order approximation as 

m 

da 

dydp T 
where 

e -E T / T =[!_(!_ q ) E T /T} 1/[1 - q) , 

e -E T / T = e -E T /T_ 

In the above equation, Ep=\Jm? +pj,, where m can be 
taken to be the pion mass m v , and we have assumed 
boost-invariance in the region near y ~ 0. The param¬ 
eter q is related physically to the power index n of the 
spectrum, and the parameter T related to mpo and the 
average transverse momentum, and the parameter A re¬ 
lated to the multiplicity (per unity rapidity) after inte¬ 
gration over pt- Given a physically determined invariant 
cross section in the log-log plot of the cross section as a 
function of the transverse hadron energy as in Fig. 1, 
the slope at large pt gives approximately the power in¬ 
dex n (and q ), the average of Et is proportional to T 
(and rriTo), and the integral over pp gives A. 

We can test the above single-particle nonextensive sta¬ 
tistical mechanical description by confronting Eq. m 
with experimental data. Fig. [8] gives the comparisons 
of the results from Eq. (ED with the experimental pp 
spectra at central rapidity obtained by different Collab¬ 
orations [65l468i |. In these calculations, the parameters of 
A, q and the corresponding n and T are given in Table 
III. The dashed line (an ordinary exponential of Et for 


1 


da 


2irpp dydpr 


= Ae- ET ' T , 


(57) 


4 We are adopting the convention of unity for both the Boltzmann 
constant kg and the speed of light c. 



FIG. 8. (Color online) Comparison of Eq. llIVTl) with the 
experimental transverse momentum distribution of hadrons 
in pp collisions at central rapidity y. The corresponding 
Boltzmann-Gibbs (purely exponential) distribution is illus¬ 
trated as the dashed curve. For a better visualization, both 
the data and the analytical curves have been divided by a 
constant factor as indicated. The ratios data/fit are shown at 
the bottom, where a roughly log-periodic behavior is observed 
on top of the g-exponential one. Data are taken from [65l - l68l . 


q —> 1) illustrates the large discrepancy if the distribu¬ 
tion is described by Boltzmann-Gibbs distribution. The 
results in Fig. 8 shows that Eq. ED adequately describes 
the hadron pp spectra at central rapidity in high-energy 
pp collisions. We verify that q increases slightly with the 
beam energy, but, for the present energies, remains al¬ 
ways q — 1.1, corresponding to a power index n in the 
range of 6-10 that decreases as a function of -y/s. 


TABLE IV. Parameters used to obtain fits presented in Fig. 
[8] The values of A is in units of GeV~ 2 /c 3 . 


Collaboration 

V'i 

A 

g 

n=l/(q~l) 

T(GeV) 

CMS (65J 

pp at 7TeV 

16.2 

1.151 

6.60 

0.147 

ATLAS [66] 

pp at 7TeV 

17.3 

1.148 

6.73 

0.150 

CMS [65] 

pp at 0.9TeV 

15.8 

1.130 

7.65 

0.128 

ATLAS [66] 

pp at 0.9TeV 

13.6 

1.124 

8.09 

0.140 

ALICE [67] 

pp at 0.9TeV 

9.95 

1.119 

8.37 

0.150 

UA1 [68] 

pp at 0.9TeV 

13.1 

1.109 

9.21 

0.154 


What interestingly emerges from the analysis of the 
data in high-energy pp collisions is that the good agree¬ 
ment of the present phenomenological fit extends to 
the whole pp region (or at least for pp greater than 
0.2 GeV/c, where reliable experimental data are avail¬ 
able) [jl|. This is being achieved with a single nonexten¬ 
sive statistical mechanical distribution with only three 
degrees of freedom with data-to-fit ratios oscillating 
about unity as in Fig. 8. Such an agreement suggests 
that the nonextensive statistical mechanical distribution 
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may not only be the phenomenological description of the 
end product of the parton showering evolution from jet 
to hadrons but may have deeper theoretical significance. 


IX. SUMMARY AND DISCUSSIONS 

Transverse momentum distribution of jets and hadrons 
provide complementary and useful pieces of information 
on the collision mechanisms and their evolution dynam¬ 
ics. The spectra of jets reveal the simple hard-scattering 
production mechanism and they carry the distinct signa¬ 
ture with a power index of n ~ 4 - 5. On the other hand, 
the spectra of hadrons contain additional subsequent dy¬ 
namics on the evolution of jets into hadrons but retain 
the power-law feature of the hard-scattering process. Re¬ 
cent description of the hadron transverse spectrum by 
a single nonextensive statistical mechanical distribution 
leads to the suggestion of the possible dominance of the 
hard-scattering process, not only in the high px region, 
but also over essentially the whole px region, for pp and 
pp collisions. The suggestion represents a synthesizing 
description linking the simplicity of the whole hadron 
spetrum for pp collisions with the production of minijets 
[3 at px of a few GeV and the production of minijets 
at low -px jUJllI into a single simplifying observation on 
the dominance of the hard-scattering over the whole px 
region in pp collisions, with a special emphasis on the 
production mechanism. 

We have searched for direct supporting evidences for 
the dominance of the hard-scattering process in the whole 
Px region at central rapidity. The first piece of evidence 
has been found by studying the power index for jet pro¬ 
duction in the lower-p^ region in the UA1 and ATLAS 
data in high-energy pp and pp collisions, where the power 
index is indeed close to 4 - 5, the signature of pQCD jet 
production. The dominance of the hard-scattering pro¬ 
cess for the production of low-px hadron in the central 
rapidity region is further supported by two-particle cor¬ 
relation data where associated particles are correlated on 
the near-side at (A <f>, Ar/j^O, with a minimum-pr-biased 
or a high-pT trigger, indicating the production of angular 
clusters in essentially the whole range of px ■ Additional 
evidence has been provided by the two-particle correla¬ 
tion on the away-side at A<p ~ tt, with a minimum-pT- 
biased or a high-p^ trigger, where a produced hadron 
has been found to correlate with a “ridge” of particles 
along At] [ 35 I. l38l - [4l| . The A cf> ~ n correlation indicates 
that the correlated pair is related by a collision, and the 
Arj correlation in the shape of a ridge indicates that the 
two particles are partons from the two nucleons and they 
carry different fractions of the longitudinal momenta of 
their parents, leading to the ridge of Ap at A <f> ~ 7 r. 

Hadron production in high-energy pp and pp colli¬ 
sions are complex processes. They can be viewed from 
two different and complementary perspectives. On the 
one hand, there is the successful microscopic descrip¬ 
tion involving perturbative QCD and nonperturbative 


hadronization at the parton level where one describes the 
detailed mechanisms of parton-parton hard scattering, 
parton structure function, parton fragmentation, parton 
showering, the running coupling constant and other QCD 
processes @]. On the other hand, from the viewpoint 
of statistical mechanics, the single-particle distribution 
may be cast into a form that exhibit all the essential fea¬ 
tures of the process with only three degrees of freedom 
000. The final result of the process may be sum¬ 
marized, in the lowest-order approximation, by a power 
index n which can be represented by a nonextensivity pa¬ 
rameter q=(n + 1 )/n, the average transverse momentum 
m,xo which can be represented by an effective tempera¬ 
ture T=mxo/n, and a constant A that is related to the 
multiplicity per unit rapidity when integrated over px- 
We have successfully confronted such a phenomenologi¬ 
cal nonextensive statistical mechanical description with 
experimental data. We emphasize also that, in all cases , 
the temperature turns out to be close to the mass of the 
pion. 


What we may extract from the behavior of the experi¬ 
mental data is that scenario proposed in [l8|, [u| appears 
to be essentially correct excepting for the fact that we 
are not facing thermal equilibrium but a different type 
of stationary state, typical of violation of ergodicity (for 
a discussion of the kinetic and effective temperatures see 
[H,[l3|; a very general discussion of the notion of temper¬ 
ature on nonextensive environments can be found in [7l|). 
It should be realized however that the connection be¬ 
tween the power law and the nonextensive statistical me¬ 
chanical description we have presented constitutes only a 
plausible mathematical outline and an approximate road¬ 
map. It will be of interest in future work to investigate 
more rigorously the stochastic parton showering process 
from a purely statistical mechanical viewpoint to see how 
it can indeed lead to a nonextensive statistical distribu¬ 
tion by deductive, physical, and statistical principles so 
that the underlying nonextensive parameters can be de¬ 
termined from basic physical quantities of the collision 
process. 


We can discuss the usefulness of our particle produc¬ 
tion results in pp collisions in relation to particle produc¬ 
tion in AA collisions. In the lowest approximation with 
no initial-state and final-state interactions, an AA colli¬ 
sions at a certain centrality b can be considered as a col¬ 
lection of binary Noi n (b) number of pp collisions. These 
binary collisions lead first to the production of primary 
particles. Successive secondary and tertiary collisions be¬ 
tween primary particles lead to additional contributions 
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in a series: 


E p dN AA /L x Ar E p dN pp A 
—-(&,P) = N hin (b)^~ - (p) 


dp 


^bin 


^inl 


(b) 


dp t dp 


dp 

dN pp dN pp E p dN(p 1 p 2 pA') 


dp { dp 2 dp 


2 

dp 3 dp 2 dp 

dN pp dN pp dN pp 


dp i dp 2 dp 3 


E p dN(p 1 p 2 p 3 -> pX') 
dp 


-m 


where E p dN(p 1 p 2 ... —i► pX')/dp is the particle distri¬ 
bution of p after binary collisions of primary particles 
Pi,p 2 , .... In addition to the primary products of a sin¬ 
gle relativistic hard-scattering EdN pp /dp represented by 
the first term on the right-hand side, the spectrum in 
AA collisions contains contributions from secondary and 
tertiary products represented by the second and third 
terms. In the next level of approximation, additional 
initial-state and final-state interactions will lead to fur¬ 
ther modifications of the ratio Rxx=dN AA /[N\, ln dN pp ] 
as a function of b and p T . 

The usefulness of our analysis arises from a better un¬ 
derstanding of the plausible reasons why the products 
from the primary pp scattering can be simply represented 
by a single nonextensive statistical mechanical distribu¬ 
tion m- For peripheral collisions, the first term of Eq. 
m suffices and the spectrum of AA collision, normal¬ 
ized per binary collision, would be very similar to that of 
the pp collision, as is indeed the case in Fig. 1 of (72|. As 
the number of binary collisions increases in more central 
collisions, the second term becomes important and shows 
up as an additional component of nonextensive statistical 
mechanical distribution with a new set of n and T pa¬ 
rameters in the region of low pp, as discussed in [73I. I74l| . 


As a concluding remark, we note that the data/fit plot 
in the bottom part of Fig. [8] exhibit an intriguing rough 
log-periodicity oscillations, which suggest corrections 
to the lowest-order approximation of Eq. m and 
some hierarchical fine-structure in the quark-gluon 
system where hadrons are generated. This behavior 
is possibly an indication of some kind of fractality 
in the system. Indeed, the concept of self-similarity, 
one of the landmarks of fractal structures, has been 
used by Hagedorn in his definition of fireball, as was 
previously pointed out in [2l| and found in analysis of 
jets produced in pp collisions at LHC This small 

oscillations have already been preliminary discussed in 
Section 8 and in [zllzzf , where the authors were able to 
mathematically accommodate these observed oscillations 
essentially allowing the index q in the very same Eq. m 
to be a complex numbeiQ (see also Refs. [78|, Q; more 
details on this phenomenon, including also discussion of 
its presence in recent AA data, can be found in |80|). 
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